Galileon gravity and its relevance to late time cosmic acceleration 
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We consider the covariant galileon gravity taking into account the third order and fourth order 
scalar field Lagrangians Lz{it) and La{tt) consisting of three and four 7r's with four and five deriva- 
tives acting on them respectively. The background dynamical equations are set up for the system 
under consideration and the stability of the self accelerating solution is demonstrated in general set- 
ting. We extended this study to the general case of the fifth order theory. For spherically symmetric 
static background, we spell out conditions for suppression of fifth force effects mediated by the 
galileon field 7r. We study field perturbations in the fixed background and investigate conditions for 
their causal propagation. We also briefly discuss metric fluctuations and derive evolution equation 
£SJ ■ for matter perturbations in galileon gravity. 

a 

■ I. INTRODUCTION 

£Nl ' The phenomenon of late time cosmic acceleration [H-O is as challenging theoretically as was the problem of black 
body radiation whose resolution unveiled many secrets of micro physics. At present, there is no definite clue for 
^ ' the theoretical understanding of the nature of cosmic repulsion. In recent years, a variety of approaches have been 
Q" 1 employed to attack the problem. According to the standard lore, the late time acceleration can be accounted for 
by supplementing the energy momentum tensor by an exotic fluid component with large negative pressure dubbed 
dark energy [5j, [6[ The simplest candidate of dark energy is provided by cosmological constant A. However, its 
- 1 . small numerical value leads to fine tuning problem and we do not understand why it becomes important today a la 
coincidence problem. 

Scalar fields provide an interesting alternative to cosmological constant though they do not address the cosmological 
constant problem. To this effect, cosmological dynamics of a variety of scalar fields has been investigated in the 
literature(see review [5] for details). They can mimic cosmological constant like behavior at late times and can 
provide a viable cosmological dynamics at early epochs. Scalar field models with generic features are capable of 
alleviating the fine tuning and coincidence problems. As for the observation, at present, it is absolutely consistent 
with A but at the same time, a large number of scalar field models are also permitted. Future data should allow to 
f^*) ' narrow down the class of permissible models of dark energy. 

It is quite possible that there is no dark energy and the late cosmic acceleration is an artifact of infrared modification 
of gravity. We know that gravity is modified at short distance and there is no guarantee that it would not suffer any 
^ , correction at large scales where it is never verified directly. Large scale modifications might arise from extra dimensional 
effects or can be inspired by fundamental theories. They can also be motivated by phenomenological considerations 
such as f(R) theories of gravity [7[ or the massive theories of gravity. However, any large scale modification of gravity 
should reconcile with local physics constraints and should have potential of being distinguished from cosmological 
constant. 

The infrared modified theories of gravity essentially contain additional degrees of freedom. The f(R) theories contain 
a scalar field which mediates fifth force and might contradicts the local gravity constraints such as the solar system 
or laboratory tests. Broadly, two mechanisms for hiding the scalar field effects locally have been employed in the 
literature. In f(R) theories of gravity, the scalar field is screened via the so-called chameleon mechanism Q, by 
making scalar field mass dependent on the local matter density. In generic models of f(R) gravity @, the chameleon 
mechanism allows to satisfy the local gravity constraints but at the same time make these models vulnerable to 
curvature singularity whose resolution requires the fine tuning worse than the one encountered in AC DM model. 
The problem can be alleviated by invoking R 2 correction but the scenario becomes problematic if extended to early 
universe @. 

An alternative possibility of large scale modification of gravity is provided by an effective scalar field w dubbed 
galileon ■ In particular such a field appears in the decoupling limit of DGP. The Lagrangian of the field respects 
the so called shift symmetry in a Minkowskian background: n — > ir + c and d^n — ¥ O^tt + where c and are 
constants. Thank to this symmetry, the equations of motion for the field contain only second derivatives. In four space 
time dimensions, there exist five Lagrangians Lj,i = 1,5 where L\ is linear in 7r, L2 contains normal kinetic term. L3 
involves three 7r's and four derivatives acting on them. This Lagrangian is obtained in the decoupling limit of DGP. 
The fourth and the fifth order Lagrangians involve four 7r's and six derivatives, five 7r's and seven derivatives acting 
on the field respectively. A general covariant form of galileon Lagrangian is obtained in Ref. (see also Ref. [3 on 
the related theme). 
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In DGP or its 4 dimensional gencralizations-galileon gravity, the effects of extra degree are suppressed using the 
Vainshtein (l3j mechanism which allows us to recover general relativity at small scales due to non-linear interaction. 
From this point of view, the DGP model is attractive model which has a self-accelerating solution, an asymptoti cally 
de Sitter solution even in the absence of vacuum energy. Unfortunately this solution suffers from instabilities [14l - fl9l | . 

galileon gravity can give rise to late time acceleration and is interesting for the following reasons: (i) It is free 
from negative energy instabilities, (ii) Unlike f(R) theories, galileon modified gravity does not suffer from curvature 
singularity, (hi) The chameleon mechanism in f(R) might come into conflict with the equivalence principle if the test 
bodies are considered as extended whereas the Vainshtein mechanism is free from this problem |20j. 

In this paper we study 4th order galileon gravity including L 3 and L4 terms in the Lagrangian. We set up FRW 
background dynamics and examine the self accelerating solution. We carry out detailed investigations on the stability 
of the solutions and discuss the spherical symmetric solutions to check the local suppression of ir effects. We also 
investigate matter perturbations in the model under consideration. 

II. LOWEST ORDER GALILEON GRAVITY AND ITS SELF ACCELERATING FRW BACKGROUND 

Recently, an interesting generalization of the DGP action in 4D was proposed in Ref.[l(|. The authors considered 
a consistent general action with a self interacting scalar field (7r) coupled. It is remarkable that the action can be 
motivated by higher dimensional considerations [2lj. In what follows we shall consider the action is invariant under 
Galilean transformation 

7r(x) -> 7r(x) + b^x 1 * + c (1) 

For the sake of simplicity, we first examine the galileon model in the lowest non-trivial order keeping up to third 
order term L 3 in the Lagrangian, 

S = J d 4 x^9 (j + cm - I (Vtt) 2 - | (Vtt) 2 ^ + S m [^ m , e 2 ^] (2) 
Similar expression occurs in the DGP model. The corresponding Einstein's equations are 

G>!/ = + cmg^v + c 2 (tt-,^-,u - ^3^( v?r ) 2 ^ + c 3 {k^k-vU-k + g^n. x n' Xp T:. p - 7r ;p [n. tli ir ;vp + 7r ;i/ 7r ;w ,]) (3) 

= PT^ + Cl + c 2 U-k + c 3 ((Dtt) 2 - TT.^TT'^ - R^lT^v) (4) 

where is the trace of the matter energy-momentum tensor, T^v = —(2/^—g) x SSm/Sg^. In spatially flat 
FRW background Eq.(j4|) gives rise to the following Friedmann equation 

3H 2 = Pm - cm + y7i 2 - 3c 3 Htt 3 (5) 
2H + 3H 2 = -cm - y 7i 2 - c 3 n 2 n (6) 

Pp m = ci - c 2 (3-Htt + + 3c 3 tt (3H 2 ir + H-h + 2Hnj (7) 
It is interesting to note that Eq.([7]) exhibits a self accelerating solution given by 

3H 2 = -cm+—n 2 -3c 3 H7r 3 (8) 

= -cm-jit 2 (9) 

which means that c\ = (we assume tt ^ 0) and H 4 = —c\/hAc\. This last condition is impossible to satisfy as c 2 
should be positive for stability of the theory. 

We therefore conclude that a stable self accelerating solution, in general, does not exist in the third order galileon 
gravity with (V7r) 2 D7r term in the Einstein frame. It is therefore necessary to invoke the higher order terms £4 and 
L 5 . In the discussion, to follow, we shall demonstrate that the desired solution can be obtained by adding the fourth 
order term in the action ^ . The analysis becomes cumbersome in the presence of 5th order term which completes 
the Lagrangian of galileon gravity. We have included the corresponding discussion and results in the Appendix. 
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III. GENERALIZATION TO NEXT HIGHER ORDER 

Let us consider the full covariant action of galileon gravity [lfj [H| • 

S = J d 4 xV=g + « + S m [^ m , e 2 ^g pv ] 

where {cj} are constants and the L? are given by 



(10) 



= 7T (11) 

L {2) = -\{Vk) 2 ^-\^ (12) 

= -i(V7r) 2 D^ (13) 

i (4) = -^(Vvr) 2 [(Dtt) 2 - n.^Tr^ + rVG^] + [U-k)^,^ - n.^n iup n> p (14) 

Varying the action fj L0[) with respect 7r and the metric g^,,, we obtained the field equation for n and Einstein 
equations 1 

c,£« = _£T (m) (15) 

<V - TW + «T$ (16) 

where f W = (l/V 3 ?) x ^ and T$ = -(2/^=5) x 6S^/Sg" u with = f d*Xi/=gL® 
where £ s and s have the following form 

= 1 (17) 

£ (2) = Dtt (18) 

£^ = (Dtt) 2 - tt^tt^ - R^^w ;v (19) 



f W = 2 (Dtt) 3 + 4 (tt./ jr./ 77./) - 6 (Dtt) (^ ;M ^ ;pl/ ) - (Dtt) (7t ;jU tt"*) i? - 2 (tt ;p tt^ tt ;i ,) fl 
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4 (Dtt) (7t ;/1 i?^ 7r ; „) + 2(Vtt) 2 (tt^ R» v ) + 8 (tt. p tt^ tt p ) + 4 tt ;i/ tt ;P(T i? pp ^) (20) 



T$ = ng^ (21) 
T$ = V,u - ^(Vtt) 2 (22) 



= 7I \m 7I > l::l7r + ft^ 7r ;A7r ;Ap 7r. p - tt p [7T ;M 7r ;l/p + 7T ; „7r ;W) ] (23) 
T$ = -4 (Dtt) tT p [tt ;p 7t ;p „ + n,„ tt ;pp ] + 2 (Dtt) 2 (7t ;p tt ;v ) - 2 (Dtt) (Vtt) 2 (tt^) - 4 (7T ;A 7T Ap tt ;p ) (*■;„„) 

+4 (tT ;A 7T ;A p,) (vT P 7T ;pl ,) - 2 (tT, Xp 7T ;Ap ) (7T ;A1 TT;,,) + 2(V7r) 2 (7T ;p ;P TV.p V ) + 4 7T :A 7T ;Ap [7r ;pp 7I> + 7I> 7T ;p ] 

+ (Dtt) 2 (Vtt) 2 ^ + 4 (Dtt) (tt-a 7T Ap tt ;p ) 5m „ - 4 (tt ; a 7T Ap 7r ;pff 7T ff ) ^ - (Vtt) 2 (Tr, prT 7T pCT ) 5pl , 
-(Vtt) 2 (tt. p *-.„) i? + i(V7T) 4 3 ^i? + 2(V7r) 2 7T p [i? PA1 tt ;!/ + R pv tt ;a1 ] - i(V7r) 4 i? p „ 

-2(Vtt) 2 (tt ;p R pa tt ;ct ) + 2(Vtt) 2 (tt p tt ct fl w „ CT ) (24) 
It may be instructive to define the effective energy density and pressure for 7r matter. Indeed, for each (i), 



1 We have T$ = -T^ 4) and £< 4 ) = -^'( 4 > compared to [H 
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VT« = n ;v £ W (25) 
which allows us to write the equation of conservation 

V^M = 7r ; „ (26) 

For each (z), assuming the perfect fluid form, we can express the field energy momentum tensor as, T$ = 
(pW + pW) + P^g^u with m m = — ex , and cr = sign(7r :0 ). The corresponding expressions for p 1 and 

P % have following form 

p W = - n p(i) = n 

p W = -i(vV) 2 PW = -i(VTr) 2 

= TrjATr^^Tr^ - (V7r) 2 D7r P^ 3 ) = T^TT^p 

p( 4 ) = 6DTTTT ;X n' Xp TT ;p - 3 (Dtt) 2 (Vtt) 2 + 3 (Vtt) 2 P< 4 ) = (Dtt) 2 (Vtt) 2 + 4D7r ;A 7r ;Ap 7r ;p - 47r ;A 7r ;Ap 7r ;p(T 7r ;,T 

+ |P (W) 4 - § (Vvr) 2 TT'PRp^ - (Vtt) 2 tt^tt^ + \R (Vtt) 4 - 2 {Vtt) 2 ^R pa ^° 

In the following section, we shall analyze the background solution of the fourth order theory. 

IV. BACKGROUND DYNAMICS 

Assuming the spatially flat background, we obtain evolution equations of the fourth order galileon cosmology, 

3P 2 = Pm + y7i 2 - 3c 3 Hn 3 + y c 4 P 2 tt 4 (27) 

2H + 3H 2 = -^ir 2 -c 3 Tr 2 n + ^c 4 n 3 (3H 2 Tr + 2Hn + 8H^j (28) 

f3p m = -c 2 (3Hn + fr) + 3c 3 tV (3# 2 tV + Hn + 2Hjtj - 18c 4 Hn 2 (3H 2 t\ + 2Hn + 3Htt) , (29) 

where we have assumed, c\ = as we do not want include the cosmological constant explicitly. In this case, the 
conservation has standard form in presence of coupling j3 

p m + 3Hp m = /3p m TT (30) 
We may also define the total energy density and pressure for the scalar field ir 

P. = y 7T 2 - 3c 3 Hn 3 + y c 4 H 2 n 4 (31) 



c 3 n 2 n - -c 4 7T 3 ^3H 2 n + 2H<h + 8H^j (32) 



which can be used to check for the total equation of state parameter io, = P^j p-^. In the next section, we discuss 
the self accelerating solution of galileon cosmology. 

V. SELF ACCELERATING SOLUTION 

A self acceleration solution is characterized by p m — and H = H = C st . 
In this case, using equation (f2"T|). we find that -k = ir = C st and 



± C 3 ± y/c 2 ~ 8C 2 C 4 

H °*° ~ 12^ (33) 

48P 2 = {*t?A± (34) 
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With A = c ^ 12c2C4±C3 V / c|-8c 2 c 4 

The existence of the self accelerating solution then implies the following conditions on constants Ci, 02,03 and c 4 , 

4~8c 2 C4 > (35) 
A+ > or A- > (36) 

It is not difficult to check the stability of the solution. Let us consider the perturbation of the two self-accelerating 
solutions, 

H = Hq + 5H, tt = 7T + Sir (37) 

It can easily be checked that SH = —3Hq5H, which means that the self-accelerating solutions are stable. 

VI. SPHERICALLY SYMMETRIC SOLUTION 

We shall now be interested in the spherically symmetric static solution . We consider a static point-like source of 
mass M, located at the origin: T^- m ' = — M5 3 (l£) and look for a spherically symmetric static solution for the field 
7r(r) described by the following differential equation 

£|JjL [rV (r)] + 2 J A [m>(rf] + 4^ A [ir'(rf] = /3M5 3 ( of) (38) 
Integration of Eq. (|3"5]) gives the following relation, 

^M^yM™)'-* ,39) 

Where r s is the Schwarzschild radius of the source. 

The conditions of existence of the solution are derived following Ref. 

if /3 > =>■ sign(c2)=sign(c4) and C3 > which means that C3 > v / 8c2~C4 if we consider the condition (|35|1 

if P < =>■ sign(c2)=sign(c4) and C3 < -\/3c2C4 which means that C3 < — v / 8c2~C4 if we consider the condition (1351) 

In case /3 < 0, at short distances, the solution is not analytic in the neighborhood of r = and we shall not consider 
this case any further. 
Whereas for j3 > 0, 



Then the galileon-mediated force is suppressed compared to the gravitational force: 



-^grav \ J \ 2/3 

At large distances, we have 



= (-) 2 <l, with rl= (^\ /2 r s (41) 



^grav ^2 

If /3 ^ C2, the galileon field can lead to late time acceleration of universe. 



Fn -2A (42) 
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VII. STABILITY 

In order to study the stability of the aforesaid static solutions, we perturb the scalar field ir: -k — > n + <fi in a fixed 
metric g pv . We have neglected the perturbations of the metric induced by the perturbations of the scalar field 0; the 
method is referred to test field approximation. 

In order to proceed with the test field approximation, let us rewrite the quadratic term in <f> in the action 

S<t> = I V^g^x jZ^-^v (43) 



with 



Zft = (44) 

zf 2 ) = -\aT (45) 

Zf" = TT^-g^Dir (46) 



J (3) 

Zfo = -2tt''' 1 R 1 ' p tt, p -2ti' u R pp tt, p ~ R^iVirf +Rtt' ix "k' v + 6Dnn^ - 6t:^ p t:.^ + 2R Pf " 7l, 7r ;p n, ry (47) 

+gT' U^p^* ~ 3(D7T) 2 + 2R pa ^ p 7T^ + ^R(Vn)A (48) 



The equation of motion for perturbations that follow from action (|43|) is 

- WZfitw - 2c l Z p \^, v + 8/3 2 0T( m ) = 0, (49) 
which we shall use in the subsequent sections. 

A. Cauchy-problem 

Following the theorem due to Leray [12], the scalar field <f> propagates causally in the effective metric G P g = —2c l Z^X 

if spacetime (M.,G£g) is globally hyperbolic. A necessary condition but not sufficient is the requirement of the 
hyperbolicity of the equation (|49| that is a Lorcntzian signature of the effective metric G P g. 

For the static spherical solution, the hyperbolicity is defined by 

c 2 + 2c 3 (27r7r + vr") + Uc^n'/r + 2ir")it r /r > (50) 
c 2 + 4c 3 7r'/r + 12c a (-k' /rf > (51) 
c 2 + 2ca(ir'/r + n") + 12c 4 7rV/r > (52) 

At large distances, we obtain the following conditions 

2 

c 2 - 36/3 2 ^>0 (53) 
c 2 r° 

c 2 + 4/?^>0 (54) 
c 2 r J 

Pa- >0, (55) 

c 2 r J 

which reduce to c 2 > at very large scales. 

At small distances, we need to impose the conditions, C4 > and C3 > 0. 
For the de Sitter phase, the hyperbolicity is defined by 
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G° c ° s = -J(^ ± +4 C2 )<0, (56) 
« 2 G C ^ = ^(A±-4c 2 ) >0, (57) 

which implies that A± > 4c 2 . 

We should however emphasize that this solution is derived when the scalar field is dominant (de Sitter phase), 
therefore any small perturbation of the scalar field leads to a perturbation of the metric and the test field approximation 
is then no longer true. 



B. Hamiltonian approach 



An alternative way to study the stability is related to the positive definiteness of Hamiltonian of the underlying 
theory. In a locally inertial frame, the Hamiltonian is 



The condition of hyperbolicity of equation (|49p is sufficient for the Hamiltonian to be bounded from below. The 
condition of hyperbolicity imposes an important restriction on sound speed which we consider next. 



C. Speed of sound 



From the equation , it is obvious to define the " sound speed" c 2 ; the condition of hyperbolicity of the equation 
restricts c s to real values c 2 > 0. It is straightforward to see that the condition of c s to be real, restrict the signature 
of the effective metric to (— , +, +, +) or (+, — , — , — ). 

However, if we also impose the positivity of the Hamiltonian, we have to consider the effective metric with the 
same signature that as that of the original metric which is (—,+,+,+), in our case . This condition for non 
superluminal behavior of the scalar field <f> is expressed by c 2 < 1. 

In case of the de Sitter phase, it is trivial to see that c 2 = ^A±+tc 2 ) < ^ (because of the conditions of stability of 
the theory (c 2 > and A± > 0)). But the problem is more delicate for the spherically symmetric solution. Indeed, 
equation (|4T))) can be rewritten as 

G%(p + Gl^d 2 (j) + G 2 ^r 2 dl(j) + first derivatives of + ... = (59) 
where 8q is the angular part of the Laplacian. 

Therefore we can define the speed of radial and angular excitations as follows, 



2 _ c 2 + 4c 3 7r7r+12c 4 ^ 2 /r 2 

Cr ~~ G™ ~ c 2 + 2c 3 (2^'/r + tt") + 12c 4 (7r' 2 /r 2 + 2^'V/r) 

2 r 2 G 2 c 2 = c 2 + 2c 3 (^7r + 7r") + 12c 4 7r'y/r 

C ° G™ ~ c 2 + 2c 3 (27r'/r + tt") + 12c 4 (7r' 2 /r 2 + 2tt'V /r) ' 

which at large distances gives rise to 
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4/3^ T JL 
c| r 3 



(62) 
(63) 



whereas for small distances, we find 
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4 = 1 (64) 
C C2 ~ 7 7 65 ) 

DC4 7T 

It is clear that at large distances, we have a superluminal behavior (c 2 > 1) of the scalar field 4> f° r the static 
spherically solution, but this behavior is physically possible if the theory does not have Closed Causal Curves (CCCs) 
which leads to paradoxes [22|, [23[ . It is known that if a spacetime is stably causal, it does not possesses CCCs which 
means that a global time can be defined. This is the case if we can define a global time for the two metrics <7 M „ and 

For the static spherically symmetric solution, we will consider the Minkowsky time rj^ v V ^tV v t = —1. 
Then 

G^fVutVvt = -ca - 2c 3 (27r7r + tt") - 12c 4 {-k' /r + 2n")n'/r (66) 
Eq. ([66|) . at large distances, reduces to 

G^V^V^^+36/3 2 ^ (67) 

which is negative iff r 6 > 36/3 2 r 2 C4/ 'c 2 . 

If this condition is satisfied then the space time (M, G c Jl) is stably causal which means that no closed timelike 
curves exist. We should emphasize that this condition is satisfied if the equation (|4Uf is hyperbolic. 

VIII. METRIC PERTURBATIONS 

Let us consider the perturbed FLRW spacetime with scalar metric perturbations in the longitudinal gauge 

ds 2 = —(I + 2(/))dt 2 + a 2 (I — 2ip) dx 2 (68) 

The linear matter perturbations S m on super horizon scales satisfy the evolution equation similar to the one in 
Einstein gravity 



6 m + 2H5 m - -^PmSm = (69) 



with the modified Newtonian constant, 



2 (c 3 tt 2 + 2/3) 2 + c 4 N 4 
4c 2 - 2c 2 tt 4 - 16c 3 iJ7r - 8c 3 7f + C4D4 



Ge // = 1 + — ^4 id^TT^ c-s- 1 „.n. ( 70 ) 



where N4 and D4 are given by 



iV 4 = 14c 2 tt 4 + c 2 tt 8 - 88c 3 #tt 5 + 4c 3 /3tt 6 + 20c 3 7T 4 7f - 64iJ/3TT 3 - 24^ 2 tt 4 + 96/37T 2 7f 

-9c 2 tt 8 - 12c 3 Hn 9 - 54c 3 n 8 n + 492H 2 n 6 - 96i^7r 7 - 48Hn 5 n + 168H-k 6 + 18/3 2 tt 8 - lU^n + 2887r 4 7r 2 



+ c 4 

+ c 2 4 [18 (11H 2 - 6Hj 7r iU + 648-ff7r y ifJ (71) 
L» 4 = -12c 2 7T 4 - c^tt 8 + 80c 3 Htt 5 - 24c 3 7r 4 if + 8 (l3H 2 + 6h) tt 2 + 96Hn^ 



+ c 4 (^9c 2 7r s + 12c 3 #7r M + 54c 3 7r s 7r - 24 [17H Z + 6HJ tt° + 288i?7r D 7fj 
+ c 2 4 (l8 (QH - 11H 2 ^ n 10 - 648£/7r 9 ^ (72) 

The study of generic models of modified gravity shows that there is a characteristic signature in the growth function 
/ = w hich can allow us to distinguish these models from ACDM and other dynamical dark energy models 

within the frame work of Einstein gravity. We expect similar features in galileon gravity. We shall address this 
important issue in our future work. 
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IX. CONCLUSION 

In this paper, we have investigated galileon gravity in its general form. The model consists of an effective field 
7r Lagrangian consisting of five terms Y^i C i-L l added to Einstein-Hilbert action such that the field equation are of 
second order. In spatially flat FRW background, we set up the evolutions equations in the model and examine the 
existence and stability of self accelerating solutions. We point out that these solutions, in general (tt ^ 0), are not 
stable in the third order galileon theory. We extend the analysis to the fourth and fifth order theory. In fourth order 
theory, self accelerating solutions exist provided that c| — 8C2C4 > and A + > or A_ > 0. We show that there is 
at least one stable self-accelerating solution in this case. The analysis is cumbersome in case of 5th order theory and 
we have included the corresponding results in the appendix. The conclusions reached in fourth order galileon theory 
are shown to hold in general. In case of the spherically symmetric static solution, we find that the solution exists 
provided that C3 > \f802C4. The solution is stable and the fifth force can lead the acceleration of the universe if we 
assume (3 ~ C2 and and C4 > 0. We find as expected that the galileon force mediated by the scalar field tt is negligibly 
small at small scales, because of the non-linear terms in the Lagrangian. However, the fifth force is of the order of 
the gravitational force at large scales in case, (3 ~ C2. 

Subsequently, we investigated the stability issues associated with the spherically symmetric solution. Using the fixed 
background method, we found superluminal behavior of perturbations as was noticed in [10j. It is really interesting 
that despite the superluminal behavior, there exist static solutions which do not possess any Closed Causal Curve 
allowing to avoid paradoxes related to micro-causality and making the solution physically acceptable. The model 
has a well posed Cauchy problem and no Closed Causal Curves exists in this model even if we have a superluminal 
behavior of the perturbation of the scalar field in the static spherically symmetric situation at large distances. 

We have included brief discussion on the metric perturbations and have set up the evolution equation for linear 
matter perturbation in the galileon gravity. In our opinion, it is important to study the growth function / = ^r^fr 
which can provide a discriminating signature of galileon gravity; we defer this analysis to our future work. 
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Appendix A: The full Lagrangian of galileon gravity: Extension of the model to the fifth order term, l/ 5 ' 
We consider the term l/ 5 ) derived in [Til] 



L< 5 



--(W) 2 [(Dirf - 3 (Dtt) (tt^ tt^) + 2 (tt./ tt./ tt./) - 3 (tt ]P tt, v tt, p<7 R^ a ) 

-18 (tt, v tt^ p R P a tt/ + 3 (Dtt) (tt, v R v p tt/ + y (V^) 2 (tt, v tt' v p tt/ R 
+3 [tt, p tt^ v w ;vp tt^ x tt, x - (Dtt) (tt.^ v tt, v(i tt/] 
(□tt) 2 (tt, p tt,,) - (tt, pv (tt, p tt^ x tt/ 



then the equations (|15I16[) are modified by 



(Al) 



| (Dtt) 4 - 15 (Dtt) 2 (tt.^ tO - y (Dtt) 2 (\7tt) 2 R - y (Dir) 2 (tt.,, tt,,) 

+20 (Dtt) {it./ it./ tt./) - y (Dtt) (tt,, tt^ v tt/ R + 15 (Dtt) (Vtt) 2 (tt, up R v ») 

+30 (On) (tt, p tt^ R vp tt/ + 15 (Dtt) (*■•„ n ;v n, pa RW°) + y ^/ - 15 (n, pv n' up n-, pa r*") 

+ y (Vtt) 2 (n, vp tt^) R + y (tt, p tt^ n, vp n/ R + y (tt^ tt'"") (tt., RP^ tt/ 
+ 15 (tt ;m tt^ 7r ; „) (^ ;pCT - 15(W) 2 (tt./ fl/ tt./) - 30 (t^ tt^ 7r w i? pCT tt ;ct ) 
-15 (tt, p tt^ R vp tt^ tt/ - y (Vtt) 2 (tt, vp tt, oX R™ px ) - 15 (tt, p tt, v TT, pa n'\ R^ x ) 

+30 (tt. a ^ tt, vp tt, R^"°) + y (W) 2 (7r iV if" tt ;p ) R - y (Vtt) 2 (tt ;w K" R prT tt/ 

-y (Vtt) 2 (tt, tt, p R aX R^ x ) + y (Vtt) 2 (tt, v tt, p R\ rX R p ™ x ) 



(A2) 



T (5) 



| (Dtt) 3 TT;,) + | (&) 3 (Vtt) V - y (Dtt) 2 (Vtt) 2 (tt^) - y (Dvr) 2 tt" [tt ;w1 7T ;!/ + tt ;p , tt ;m ] 

+ y (Dtt) 2 (tt, p tt^ tt/ g pv + 15 (Dtt) (Vtt) 2 (tt, p<7 tt'/) - 15 (Dtt) (tt, p it" tt/ (tt,/ 
15 

-y (Dtt) (TT, pa TT'' pa ) (tt, p tt/ + 15 (Dtt) (tt" tt w ) (tt''' 7 tt,/ + 15 (Dtt) tt, p n' p ° [tt, gp tt, v + tt, ov tt/ 
-y (Dtt) (Vtt) 2 (tt ;ctA 7T ctA ) fl/M , - 15 (Dtt) (tt, p tt^ tt, gX tt/ g pv - y (Dtt) (V^) 2 (tt, p tt/ R 
+ y (Dtt) (Vtt) 2 ^ [R ap tt,, + R av tt/ - y (Dtt) (Vtt) 2 (tt ;ct iT A 7r ;A ) 5pt/ 
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+ y (Dtt) (VVT) 2 (tT" 7T A i?^) + y (W) 2 (lT, aX ^ A ) (tT ; ^) - 15(W) 2 (tT^ TT^ 7T ;Ai ,) 

+ 15 (tt ;p ir> pa tt ;(7 ) (tt ;mA 7T A „) + 15 (tt ;p Tr" 17 7t ;ctA tt a ) (tt ; ^) + 5 (tt./ tt./ tt./) (7T ;jU 7i>) 
15 

+ "^" (''"lO'A 7r ,<T ) 7I" P 71> + 7I";piy 7i" ;A ,] — 15 7T P 7T ;pCT 7T ,<T [7T ;Ap 7T : „ + 7T ;A „ 7T ;p ] 

15 

-15 7T ;p 7T ;pA 7T i<T [tT ;Am 7T ;(TV + 7T :A „ 7T ;(Tp ] + 5(Vtt) 2 (tT ;(T A 7T./ 7T./) ff/4w - y (tT ;p 7^ 7!>) (tT ;Ak 7T Ak ) ^ 

15 15 

+ 15 (7r ;p Tr'" 7 7T ;ffA TP*" 7T ;K ) 5ply + — (Vtt) 2 ^^ [n, afl 7T ;i , + 7T ;w 7T ;M ] R - y (Vtt) 2 (tT ;(T TT^ 7T ;A ) i? 5pl/ 

+ y(V^) 2 tt ctA tt ; a) + y (W) 2 (tt,., iT A tt ;A ) (tt ; ^) + y (Vtt) 2 (tt,^ R° x ) (tt ;p , tt ;!/ ) 

-y(V7r) 2 7r ;ff 7T ;ffA [i? Ap 7r ;l/ + i? A „ 7T ;ai ] - y (V7t) 2 7T ;A 7T ;ff TT.^ + fl A v 7T;a M ] 

-y (V7T) 2 7T ;CT i? CTA [7T ;Ap 7T ; „ + 7T ;A „ 7T ;p ] + 15(Vtt) 2 (n ;r7 TT° X R Xk 7T' K ) <? pt , 

-y (Vtt) V" 7 TT^ [i? pAffK n- v + i?, AffK 7T ;/ J + y (V7r) 2 7T ;CT 7T A [i? pCT A K TT^ + R vaXtt TT^] 

-y (VttJV tt ;ctA 7T K [i^ Ai , K + i?„ VK ] + y (V^) 2 (7T ;(7 tt ;A ^ ;kt iT KAr ) g^. (A3) 
The Friedmann equations for this model are 



3H 2 = Pm + |^ 2 - 3c 3 H^ + y^ff 2 * 4 - y^ff 3 * 5 (A4) 
2H + 3iJ 2 = ~y7T 2 - c 3 7r 2 ^ + ^c 4 tt 3 (3H 2 ir + 2ij7i + 8i27r) - yc 5 iJ7r 4 (2iJ 2 7r + 2£/^t + sW) (A5) 

/9pm = -c 2 (3i?Tr + if) + 3c 3 tt (3iJ 2 7r + #tt + 2iJif) - 18c 4 -ff7r 2 (3iJ 2 7r + 2iJ7r + 3//"tt) 

75 / \ 

+ yc 5 ij 2 i 3 (3iJ 2 7r + 3i77r + 4Hn) (A6) 

Therefore the self accelerating solution exist if there is a real solution of the equation 

75 

c 2 - 3c 3 A + 18c 4 X 2 - — c 5 X 3 = (A7) 
c 2 - 9c 4 X 2 + 30c 5 A 3 < 0, with X = i7 vr (A8) 

If we have a solution of this system therefore the self accelerating solution is stable. In fact if we consider a 
perturbation of the self accelerating solution H = Ha + SH and n — ttq + 5ii 1 it is straightforward to see that, 
SH = -3H 5H. 

We found that the spherical symmetric solution is not modified by the fifth term as it was noticed in Ref. [To[ ■ 



